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ABSTRACT 

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main 
sources have been suggested: asymptotic giant branch (AGB) stars and massive stars rotating near the break-up limit (spin stars). We 
propose massive binaries as an alternative source of processed material. We compute the evolution of a 20 M Q star in a close binary 
considering the effects of non conservative mass and angular momentum transfer and of rotation and tidal interaction to demonstrate 
the principle. We find that this system sheds about 10 M Q of material, nearly the entire envelope of the primary star. The ejecta are 
enriched in He, N, Na, and Al and depleted in C and O, similar to the abundance patterns observed in gobular cluster stars. However, 
Mg is not significantly depleted in the ejecta of this model. 

In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it 
remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of 
stars. We estimate that the amount of processed low-velocity material ejected by massive binaries is greater than the contribution of 
AGB stars and spin stars combined, assuming that the majority of massive stars in a proto-globular cluster interact with a companion 
and return their envelope to the interstellar medium. If we take the possible contribution of intermediate mass stars in binaries into 
account and assume that the ejecta are diluted with an equal amount of unprocessed material, we find that this scenario can potentially 
provide enough material to form a second generation of low-mass stars, which is as numerous as the first generation of low-mass stars, 
without the need to make commonly adopted assumptions, such as preferential loss of the first generation of stars, external pollution 
of the cluster, or an anomalous initial mass function. 
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1. Introduction 

For a long time star clusters have been considered as ideal- 
ized single-age, chemically homogeneous stellar populations. 
However, it has recently become clear that many clusters show 
multiple main sequences and sub giant bran ches and extended 
horizontal branches (e.g. iPiotto et al.l [20071) . implying the ex- 
istence of multiple populations within one cluster.In addition, 
large star-to-star abundance variations are found for light ele- 
ments such as C, N, O, Na, and Al, while the composition of 
heavier elements (Fe-group and a-elements) seems to be con- 
stant. Field stars with the same met allicity do not exhib it these 
abundance patterns (for a review see iGratton et al1l2004 . These 
chemical variations have been interpreted as originating from the 
presence of both a "normal" stellar population, exhibiting abun- 
dances similar to field stars of the same metallicity and a sec- 
ond population of stars formed out of material processed by hy- 
drogen burn ing via the CNO-cycl e and by the Ne Na and MgAl 
chains (e.g. iPrantzos et al1l2007l) . According to ICarretta et al.l 
(120091) . 50-70% of the stars in gloular clusters belong to the sec- 
ond population. 

Two sources of processed ejecta have been proposed: the 
slow winds of massive AGB stars, which enrich their connec- 
tive envelopes with H-burning products (| Ventura et all 1200 ll 
lD'AntonaetaIll2002l: iDenissenkov & Herwigl 120031) and fast- 
rotating massive stars (we refer to these as spin stars), which 



are believed to expel processe d material centrifugally when 
they reach break-up rotation dPrantzos & Charbonnell 120061; 
iDecressin et al.ll2007b|) .In this scenario a first generation of stars 
is formed out of pristine material. Their low-velocity ejecta 
are trapped inside the potential well of the cluster and provide 
the material for the formation of a second generation of stars. 
Although both proposed sources are promising, matching the ob- 
served abundance patterns and providing enough ejecta for the 
formation of a second generation that outnumbers the first gen- 
eration have proven to be two major challenges. In this Letter we 
propose massive binaries as a candidate for the internal pollution 
of globular clusters. 

2. Binaries as sources of enrichment 

Interacting binaries can shed large amounts of material pro- 
cessed by hydrogen burning into their surroundings. A clear 
example is the well-studied system RY Scuti. It is undergoing 
rapid mass transfer from a 7 M supergiant to its 30 M com- 
panion. Mass is lost from the system via the outer Lagrangian 
points into a circumbinary disk and a wider double toroidal neb- 
ula. The nebula shows signatures of CNO processing: it is en- 
riche d in helium and n i trogen and depleted in o xygen and car- 
bon dSmithet all 120021; iGrundstrom et al J 120071) . In contrast to 
the high-velocity radiatively driven winds of massive stars, these 
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Fig. 1. Composition of the slow ejecta of the modeled binary system (Sec. 3) as a function of the ejected amount of mass. The 
mass fraction X of the main stable isotope of each element is given relative to the initial mass fraction X{, except for Mg where we 
added 24 Mg, 25 Mg and 26 Mg. The average X av and the most extreme mass fraction X ex are given in each panel on a logarithmic scale: 
[X] = log l0 (X/Xi). For helium we show the absolute mass fraction Y instead. Mass ejected during the first and second mass transfer 
phases is separated by a thin vertical line. 



ejecta have low velocities. ISmith et ail d200ll) measure expan- 
sion velocities ranging from 30 to 70 km s in the nebula of 
RY Scuti, which are lower than the present-day escape velocity 
of massive globul ar clusters. Further more, the nebula shows ev- 
idence of clumps (ISmith et ai1l2002l) and dust in the outer parts 
dGehrzetal.1120011) . which may serve as seeds for forming a sec- 
ond generation of low-mass proto stars. 

Evidence of severe mass loss from interacting binaries 
comes from a wide variety of observed interacting and post- 
interaction systems. This appears to be a common phe- 
nomenon for many interacting binaries. Various authors have 
inferred highly non conservative evolution for Algols, sys- 
tems that are currently undergoing stable mass transfer (e.g. 
Refsdal et all 1 19741 Ide Greve & Linnelll 1 1 994t iFigueiredo et all 



19941: Ivan Rensbergen et all 120061: |Pe Mink et ail 120071) . Most 
notable are short-period binaries containing a compact object, 
e.g. cataclysmic variables, X-ray binaries, binary radio pulsars 
and double white dwarf systems. Their formation requires a 
phase of severe mass and angular momentum loss by ejection 
of a common envelope. Direct evidence of this type of evolution 
comes from planetary nebulae with close binary nuclei, which 
appear to have recent ly emerged from th e common-envelope 
phase (for a review see llben & Livio|[T993l) . 

Theoretical considerations support the idea that most inter- 
acting binaries shed large amounts of mass. Three-dimensional 
hydrodynamical simulations of the mass transfer stream and 
accretion disk of the interacting binary B Lyrae predict that 
50% of the transferred mass is lost (Bisikalo et al. 2000; 
iNazarenko & Glazunoval 12006b . In addition, lUlrich & Burgerl 
( 1 19761) showed that the accreting star is driven out of ther- 
mal equilibrium and expands. This can lead to contact and 
to strong mass and angu lar momentum loss from the system 
jFlannerv & Ulrichl 1 1977b . Furthermore, iPacketl dl98ll) noted 
that the accreting star reaches break-up rotation after gaining 
only a few percent of its own mas s. Rapid rotation is found 
for many accreting stars in Algols dBarai et al. 1 12004b and this 
mechanism has been proposed to explain the formation of Be-X- 



ray binaries (e.g. lPols et al.lfl99~il) . In principle, tides can coun- 
teract the effect of sp in-up by mass transfer in close binaries. 
iPetrovic et al.l d2005al) have computed detailed binary evolution 
models taking these effects into account. They find that massive 
binaries with initial periods as short as 3-6 days lose 70-80% of 
the transferred mass on average. For wider and more massive 
systems, they expect even less conservative mass transfer, such 
that nearly the entire envelope of the primary is returned to the 
interstellar medium. 

3. Composition of the ejecta 

To investigate the yields of a typical massive binary, we 
employ a stat e -of-the- art bi nary evolut i on co de described by 
IPetrovic et ail d2005bl) and lYoon et al l (120061) . The effect of 
mass and angular momentum loss on the binar y orbit is com- 
puted according to [Podsiadlowski et ail d 1992b . with the spe- 
cific angular momen t um of the wind calculated according to 
iBrookshaw & Tavanil d!993b . Mass transfer is modeled follow- 
ing iRitterl (Il988l). Ti dal interaction is modeled as described in 
iDetmers et al.l ( 2008b . Non conservative mass transfer is mod- 
eled self-consistently: it results from the interplay between spin- 
up by mass transfer, tidal interaction, and rotationally enhanced 
mass loss. To follow the nucleosynthesis up to the advanced 
stages of hydrogen burni ng, we updated the r eaction rates to the 
NACRE 99 compilation (lAngulo et al.lll999b . As initial compo- 
sition we assumed an g- enhanced mixture with a metallicity of 
Z = 5 x 10~ 4 , following Ibecressin et all d2007bb . We assumed 
masses of 20 and 15 M for the two stars, an orbital period of 12 
days and initial rotation rates synchronized with the orbital rev- 
olution. Because of these low initial rotation rates, the effect of 
rotationally induced mixing on the low-velocity ejecta is negligi- 
ble. We follow the evolution from the onset of hydrogen burning 
until central carbon burning. 

Shortly after hydrogen exhaustion in the center, the primary 
expands and starts to transfer mass to its companion. Initially 
the secondary star efficiently accretes all the transferred mass 
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and associated angular momentum and spins up. After accreting 
about 1 .5 M it approaches critical rotation. From this moment 
on, the majority of the transferred mass is ejected from the sys- 
tem. 

The transition from conservative to non conservative mass 
transfer, when the accreting star reaches critical rotation, pro- 
vides an interesting selection mechanism. Initially, when the out- 
ermost unprocessed layers are transferred, the companion star 
efficiently accretes all the material. By the time deeper layers of 
the donor star, which do show signatures of nuclear processing, 
are exposed, mass is ejected from the system. 

After transferring nearly its entire envelope, the donor con- 
tracts, ignites helium, and becomes a Wolf-Rayet star. In the 
meantime the secondary spins down because of angular momen- 
tum loss in its rotationally enhanced wind. A second phase of 
mass transfer sets in when the Wolf-Rayet star expands during 
He-shell burning. This time about 1 M Q is transferred, which is 
initially accreted by the secondary but ejected shortly afterwards 
by its rotationally enhanced wind. Our computation ends after 
the primary star ignites carbon and fills its Roche lobe a third 
time. Soon after, it will explode as a Type Ib/c supernova. 

Figure Q] shows the composition of the material ejected from 
the system during the two phases of mass transfer. The first 
2 M are relatively unprocessed and resemble the pristine com- 
position, except for the depletion of fragile elements such as 
lithium (not plotted). The next 2 M of ejected material are pro- 
cessed by CN cycling (nitrogen being enhanced up to a fac- 
tor five), followed by about 4 M showing He enrichment and 
the O-Na anti-correlation. After the ejection of 8 M , a sudden 
change in slope is visible for all elements except carbon. The 
layers of the donor star, which are exposed at this moment, were 
part of the convection zone above the H burning shell. Here, the 
temperatures were high enough for proton captures onto 25 Mg 
and 26 Mg, leading to an increase in the aluminum abundance 
by a factor two. Between the first and second mass transfer 
phases, the primary star loses mass in the form of a fast Wolf- 
Rayet wind, and even deeper layers of the primary are exposed. 
Because the these high-velocity winds are likely to escape from 
the cluster, we have exclude the mass ejected during this phase 
in Fig.Q] 

The average composition of the ejecta are c omparable to the 
yields o f a 60 M spin star model computed by iDecressin et alj 
d2007bl) by adopting the same set of reaction rates. For example, 
we find an average helium mass fraction of 0.30 and an enhance- 
ment in sodium by 1 .0 dex compared to 0.32 and 1.3 dex, respec- 
tively, in the 60 M spin star. The sum of carbon, nitrogen, and 
oxygen is constant within a few percent, consist ent with obser- 
vations (see references in IDecressin et al.ll2.009l) . However, the 
temperatures are not high enough for efficient proton captures 
onto the most abundant isotope of magnesium, 24 Mg, adopting 
the recommended reaction rates. The onset of this reaction is 
required to explain the full range of aluminum abundances ob- 
served in some globular clusters. We expect that considering 
more massive or wider binaries (in which mass transfer starts 
in a later evolution stage) might alleviate this problem. 

In shorter-period binaries, tides can counteract the effect 
of spin-up: more mass is accreted before the accreting star 
reaches break-up rotation and less mass is lost from the system. 
However, the ejecta originate from the last exposed (deepest) 
layers of the donor star and will therefore show more strongly 
pronounced anti-correlations. For initially wider systems, we ex- 
pect that nearly the complete envelope is ejected. In addition 
stars may have considerable rotation rates at birth. For these 
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Fig. 2. Mass-weighted lKroupal d2001l) IMF as a function of stel- 
lar mass. The surface areas indicate the mass contained in the 
first generation of long-lived low-mass stars (dark blue), the 
ejecta of AGB stars (dotted line), spin stars, i.e. fast-rotating 
massive stars (dashed line), and massive (red) and intermediate- 
mass (orange) binaries. Percentages indicate the fraction of mass 
relative to the total mass contained in stars of the first generation. 
See Sect. |4]for details. 

stars, mixing induce d by rotation can lead to processing of t he 
whole envelope (e.g. lYoon et al.|[2006i : IDecressin et al.ll2007bl) . 

Besides the effects discussed above, the (remainders of the) 
binary can still shed large amounts of H-processed material 
at low velocities in various ways. The companion star is now 
rapidly rotating and may reach break-up rotation again towards 
the end of its evolution as proposed in the spin star scenario. 
Processed mass can be ejected during a phase of reverse mass 
transfer from the secondary onto the compact object, if the sys- 
tem remains bound after the supernova explosion of the primary 
star. 

4. Mass budget 

One of the main challenges for the two previously proposed 
sources of pollution, massive AGB and spin stars, is to provide 
the large amount of ejecta needed to create a second popula- 
tion that is larger than the first population. The population of 
low-mass stars (0.1-0.8 M ), which can still be observed today, 
represents 38% of th e stellar mass i nitially present in the cluster 
assuming a standard [Kroupal d2001l) initial mass function (IMF) 
between 0.1-120 M , see Fig [2] The ejecta of AGB stars with 
initial masses between 4 and 9 M represent up to 8.9% of the 
initial stellar ma ss (assuming an initial-final mass relation by 
ICiotti et alii 19911) . For spin stars this fraction is 3.4%, if one as- 
sumes that every massive star is single and born with a rotational 
ve locity high enough to r each break-up rotation (using models 
bv lDecressin et al.ll2007bl) . There are not enough of these ejecta 
to create a second generation which is equally numerous as the 
first generation, even when we assume that the second genera- 
tion consists only of low-mass stars and that star formation is 
very efficient, see Fig [2] 

Two rather extreme solutions have been proposed. (1) The 
IMF was highly anomalous, favoring the formation of the pol- 
luting stars with respect to the long-lived low-mass stars that we 
observe today. Even tho ugh we have no direct constraints on the 
IMF of globular clusters. lKroupal d2002l) finds that the IMF is re- 
markably uniform in stellar populations with very different prop- 
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erties. (2) Clusters were initially at least 10-20 times more mas- 
sive and they have primarily lost low-mass stars from the first 
generation as a res ult of the dynamical evolution and tidal strip- 
ping o f the cluster dDecressin et al.ll2007al 120081: iD'Ercole etail 
120081) . In this section we investigate to what extent the ejecta of 
massive binaries can alleviate this conundrum. 

Even though the current f raction of detected binaries in glob- 
ular clusters is not high (e.g. lDavis et alj|2008l) . this is not nec- 
essarily thecasefor thehigh-mas^_s^^ in the 
cluster. ISana et al.l d2008l) and iMason et all d2009t) find a min- 
imum binary fraction of 60-75% for O stars associated with 
clusters or OB associations. In globular clusters, these frac- 
tions may even be higher. In this environment, close binaries 
can be created and tightened during and after the star forma- 
tion pro cess, for example by the dissipative interaction with 
gas (e.g. lBonnell & Batg|2005l) and by three-body interactions, 
such as the Kozai mechanism in combination with tidal friction 
dFabrvckv & Tre maine 2007). With massive stars preferentially 
residing in the dense core of the cluster, where the dynamical en- 
counters are most frequent, it is not unreasonable to assume that 
the large majority of massive stars interact by mass exchange. 

Let us assume that every massive star is a member of an inter- 
acting binary. In Sects. 2 and 3, we argued that nearly the entire 
envelope of the donor is returned to the ISM. For simplicity, we 
neglect the contribution of the secondary star after it has been 
spun up by mass transfer or during a possible phase of reverse 
mass transfer, and we assume that the entire envelope of the pri- 
mary becomes av ailable for star formation. We a ssume helium- 
core masses as in iPrantzos & Charbonnell (12006b for stars more 
massive than 10 M . Under these assumptions, the slow ejecta 
of massive binaries represent 13% of the mass originally present 
in stars: more than the ejecta of AGB and spin stars combined. 

Measurements of lithium suggest that the ejecta of the firs t 
generation are diluted with pristine gas dPasquini et al.l 120051) . 
Together with an equal amount of pristine gas, the ejecta of bi- 
naries with donors more massive than 10 M would represent 
26% of the initial cluster mass (compared to 38% contained in 
the first generation of low-mass stars). The adopted lower mass 
limit for our binary scenario is rather arbitrary. If we take the 
potential contribution of intermediate mass stars (4-10 M ) into 
account according to this scenario, there would be enough ejecta 
to form a second population of chemically peculiar stars that out- 
numbers the first generation of normal stars. The assumptions in 
this scenario can be relaxed if the evaporation of stars from the 
cl uster primarily affects the fi rst stellar generation, a s suggested 
bv lD'Ercole et alj (120081) and lDecressin et al.l d2008l) . 



5. Conclusions 

We propose massive binaries as a source for the internal pollu- 
tion of globular clusters. The majority of massive stars are ex- 
pected to be members of interacting binary systems. These re- 
turn most of the envelope of their primary star to the interstellar 
medium during non conservative mass transfer. We show that 
there may be more polluted material ejected by binaries than by 
the two previously suggested sources: massive AGB stars and 
the slow winds of fast-rotating massive stars. After dilution with 
pristine material, as lithium observations suggest, binaries could 
return enough material to form a chemically enriched second 
generation that is as numerous as the first generation of low- 
mass stars, without the need to assume a highly anomalous IMF, 
external pollution of the cluster or a significant loss of stars from 
the unenriched first generation. 



In addition to providing a new source of slowly-ejected en- 
riched material, binary interaction also affects the previously 
proposed scenarios. Binary mass transfer naturally produces a 
large number of fast-rotating massive stars that may enrich their 
surroundings even more. Binary interaction will also affect the 
yields of intermediate-mass stars. Premature ejection of the en- 
velope in 4-9 M stars will result in ejecta with less pronounced 
anti-correlations, as suggested in the AGB scenario. On the other 
hand, we expect that binary-induced mass loss may also prevent 
the dredge-up of helium-burning products. 

For a detailed comparison of the chemical predictions of 
this scenario, binary models for a range of masses and orbital 
periods are needed and population synthesis models are essen- 
tial to fullly evaluate the mass budget of the different sources. 
Finally, some peculiar feature, such as the apparent presence of 
distinct, chemic ally homogene ous subpopulations in a> Cen and 
NGC 2808 (e.g. lRenzinill2008l) deserves further attention. 
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